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Abstract 

It is shown that colour transparency causes nonperturbative colour-singlet final-state interactions to 
have a negligible effect on the production rate and the dijet mass spectra in e + e~ — > WW. However, 
the same cannot be said of nonperturbative colour-octet exchange, for which we show that there are 
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£T) ■ indications of observable effects, though we are unable to present precise estimates. 
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O 1 Introduction 

An important aim at LEP2 is to measure the mass of the M^-boson to high accuracy, perhaps to within 
50 MeVW. A preferred method, because it offers the prospect of the most statistics, is to produce a 
pair of W's, each of which decays into a pair of quark jets. The invariant mass of each jet pair is then 
measured. 
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An obvious question is whether final-state interactions among the jets will cause a problem: the W- 
particles have a short lifetime, so the quarks are close together when they are produced and interactions 
among them could well be significant. There have been several calculations of colour exchange among 
the quarks M , mostly reaching the conclusion that the effect is quite small. 

In this paper, we first examine the effect of colour-singlet exchange. Because the relative energies 
of the quarks are large, we use soft-pomeron phenomenology M to model this exchange. While this 
phenomenology is well-established and has had many successes, its extension to the present problem 
does involve uncertainty. Nevertheless, our conclusion that there is a negligibly small effect is prob- 
ably reliable. This result comes about because colour-transparency effects^] suppress colour-singlet 
exchange. 

We then go on to perform a similar calculation of colour-octet exchange. We model this exchange 
using Cornwall's solution to the Schwinger-Dyson equations^], which provides a well- motivated way 
to handle the nonperturbative region of gluon exchange. In the case of perturbative gluon exchange, 
there is an infrared divergence, corresponding to soft exchange, which is cancelled by adding in the 
contribution from soft gluon emission. However, these infrared divergences are not genuine, because 
the nonperturbative corrections to the propagator, which have their origin in confinement, remove 
them. The Cornwall formalism does not provide a method of calculating the gluon emission, but since 
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Figure 1: e + e~ — > 4 quark jets via WW, (a) the Born diagram and (b) with final-state interaction. 

it effectively gives the gluon a mass, any emission that does occur cannot be soft. Hence there is 
nothing to cancel the correction to the cross section from soft gluon exchange. We therefore calculate 
soft nonperturbative gluon exchange and maintain that it is likely to provide a lower bound to the 
true correction to the cross section from colour-octet exchange. We find that this lower bound is far 
from small. 



Formalism 

We consider the differential cross-section d 2 a/dM\ dM 2 for the e + e~ annihilation into WW and the 
decay of each W into a quark-antiquark pair, WW — ► qqqq, with Mi and M 2 the invariant masses 
of the jet pairs. One of the Born diagrams for this process is shown in figure la. We concentrate on 
the corrections to this diagram; we do not expect any significant difference in the result for the other 
Born diagram, in which the e + e~ transform into WW through neutrino exchange. 

We suppose that the interactions among the final quarks are two-body interactions, as depicted in 
figure lb. We suppose further that they do not flip the helicities of the quarks. Then the ratio of the 
amplitude of figure lb to that of figure la is independent of the spin states of the leptons and the 
quarks and is 

R(P p „ „ ~ M$y)(P$ ~ M$y) T qq ( Pl , P2 ,A) 

R(P 1 ,P 2 , Pl ,p 2 )- l — A Jd A ((Pi + A)2 _ M ^ )((p2 _ A)2 _ M ^ (1) 

Here the W mass is complex, My/ = M — iT. The quark-quark scattering amplitude T qq includes the 
two off-shell lines p\ = p\ + A and p 2 = p 2 — A. It is averaged over the quark spin states. We have 
to sum (1) over the four different ways of attaching the interaction to the quark and antiquark lines. 

We define v = P\ ■ P 2 . It is useful to introduce linear combinations Ri and R 2 of Pi and P 2 satisfying 
R\ = = R 2 '. 

P 1 =R 1 + \ X R 2 p 2 = R 2 + \ 2 R X (2a) 



with 



V - yV - M\Ml _V-y/lS 2 - MjMl 



Then vr = R\- R 2 satisfies vr = u(l + A1A2). We parametrise 



a 

A = —R 1 -^—R 2 + 5 (3a) 
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where <5 • R\ = = 8 ■ P 2 , so that S is a two-dimensional anti-euclidean vector. Then 

J d 4 A = j dudf3d 2 5 (36) 

We find that it is consistent to assume that most of the contribution to the A integration will arise 
from values of a and (3 that are much less than ur, so that 

t = A 2 ~ 5 2 (4a) 

and the squared 4- momenta of the W's are 

P 2 = (Pi +A) 2 ) ~ A x -/3 + Aiq 

P| = (P 2 - A 2 ) ~ A 2 - a + A 2 /3 (46) 

with 

Ai = Ml + 5 2 A 2 = Mf + 5 2 (4c) 

With 

pi = xRi + y'P 2 + tti P2 = x'Ri + yP 2 - 7r 2 (5a) 

where 7Ti and 7r 2 are each transverse to both Pi and P 2 and so again are two-dimensional anti-euclidean 
vectors, 

j d 4 Pl 5 + {p\) o+((Pi - pi) 2 ) = ^- J dxdyi d 2 7T, S(y' - X 1 (l - x)) 5(tt 2 + \ lX (l - x)) 

J d 4 P2 5 + (p 2 )7r + ((P 2 -p 2 ) 2 ) = ^-J dx'dyd 2 7T 2 5(x' -X 2 (l-y))5(7r 2 + \ 2 y(l-y)) (56) 
Further, the energy variable for the qq interaction is 

v = Pi -P2 = v R (xy + AiA 2 (l - x)(l - y)) -t\ x --k 2 (6a) 
and the squared 4-momenta of the virtual quarks are 

pi = (pi + A) 2 ~ Pi - px + Ai«(l - x) 

p 2 = (p 2 - A) 2 ~ B 2 - ay + A 2 /3(l - y) (66) 

with 

Pi = 5 2 - 2tti • 5 B 2 = 5 2 - 2tt 2 ■ 5 (6c) 

The og-interaction amplitude T is a function of the energy z>, the momentum transfer t, and the 
squared 4-momenta p 2 and p 2 of the virtual quarks. According to standard analyticity properties^], 
the singularities of T are confined to the upper halves of the complex planes of the two variables (6b). 
This is conveniently expressed by the representation 

/■OO 

T= dKxdK 2 T(v,i, K\, k 2 ) expi(K\p 2 + K 2 p 2 ) (7a) 
Jo 

We introduce the W propagators 
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Figure 2: Regions of the (x, y)-plane for which the integral (8a) is nonzero. The curve is E = 0, with 
E given in (8c) 

insert the expressions (4b) and (6b) into (7) and perform the integrations over a and (3 by closing 
each integration contour in the appropriate half of the complex plane: 

dad(3 TV1V2 = - / dK 1 dK2 / T(i>,i,K 1 ,K2)- r— — 0(/«i + Gk 2 ) 6(Fk 1 + gn 2 ) expz(uiKi +u 2 k 2 ) 

Jo 1 - AiA 2 

(8a) 



where 



u 1 = B 1 + (Ai - M$y)f + (A 2 - M^)F 

/ = A 1 A 2 (1-Ax)/(1-A 1 A 2 ) 
F = Ai(l-2x)/(l- AiA 2 ) 

A = 



u 2 = £ 2 + (Ax - M^)G + (A 2 - M^) 5 

= A 1 A 2 (1-Ay)/(1-A 1 A 2 ) 
G = A 2 (l — 2y)/(l — AiA 2 ) 



l + AiA 2 
AiA 2 



(86) 



The two ^-functions in (8a) give different limits to the K\ and k 2 integrations for different ranges of 
values of x and y, corresponding to different signs for f,g,F,G. This is shown in figure 2, where the 
numbers denote the regions of the (x, y)-plane for which the integral is non-zero. The region below 
the dashed curve is the region E < 0, where 

E = fg-FG (8c) 

Suppose, as an example, that T(u,t,p 2 ,p 2 ) were to have the simple factorised form 

1 



p 2 — n 2 



(9) 



Then the integral (8a) is (2tt) 2 A(i>,i)/(\ — \\\ 2 ) times the following: 
Region 1: F,G > 0, f,g<0 

E9{-E) 1 
FG D 1 D 2 
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Region 2: F, G > 0, / > 0, g < (so that E < 0) 

1 
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Region 3: F, G > 0, / < 0, 5 > (so that E < 0) 
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Region 4: F,G > 0, /, 5 > (so that E < 0) 
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Region 5: F > 0, G < / > 0, c/ < 

-9(E)— ^- 9(—E)- '' 



Region 6: 7 1 < 0, G > / < 0, # > 



-"( £ »^- 9 <- £ )75k (10<1) 

where 

U 1 =u 1 - /if U 2 = u 2 - nl 

D 1 = U 1 + CU 2 
D 2 = dJ r + U 2 

C = -f/G c = -g/F (106) 

This complexity is not unexpected. In coordinate space there are many different contributions, cor- 
responding to the possible different time-orderings of the vertices in figure lb. Although we cannot 
identify a direct correspondence between the different time-orderings and the various cases in (10), we 
speculate that there is a connection between the two. 

For both singlet and octet exchange we perform the integrations over the two components of 5 nu- 
merically. In each case we find that there is an important contribution from values of x and y that lie 
near the curve E = 0. When E = 0, 

GD X = GBi - fB 2 FD 2 = FB 2 - gB x 

and so the dependence of D\ and D 2 on My/ disappears. So, on the curve E = in region 1, the 
amplitude of figure lb no longer depends on the W mass: presumably this corresponds to the Ws 
decaying very quickly, before they can propagate, so that the quarks are still close together and their 
interaction is therefore enhanced PL 



Soft-pomeron exchange 

In the case of colour-singlet exchange, we find that the amplitude is so strongly peaked near to E = 
that we do not need to consider interference between the different ways of attaching the exchange to 
the quark lines. We calculate R, defined in (1) and appropriately summed over the different ways of 
attaching the exchanges to the quarks, and integrate |1 + R\ 2 over the angle between tt and it' . So we 
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include the square of figure lb, and also the interference with the Born term of figure la. In this way 
we determine the ratio 



4 ^CORRECTION 



K(M 1 ,M 2 ,x,y) = 



4^BORN 



d 4 <7 



(11) 



dM 1 dM 2 dxdy I dM 1 dM 2 dxdy 
where x and y are defined in (5a). 

We now introduce the specific soft-pomeron-exchange form for the final-state interaction. ThenW 

A{v,i) = 2Pl{2vT^ 
a(t) = 1 + e + a'i 

l3 = 2 GeV -1 e = 0.08 a' = 0.25 GeV -2 (12a) 

We need a form factor for the coupling of the pomeron to the off-shell quarks. The best available 
choice to use (9) and replace the function 4>(p 2 ,fi 2 ) with 



0(^,0) -<f>(p 2 ,ri) 



(126) 



with fiQ lGeV. 



The result of our numerical computation of 1Z, defined in (11), is that it is extremely small over 
almost the whole (x,y) plane, at the one per mil level or less. This is to be attributed to colour 
transparency^ , which in our calculation manifests itself as a strong cancellation between the two 
terms in (12b). 



Octet exchange 

According to Cornwall's calculation^] , colour-octet exchange between quarks can be well approximated 



by 
with 



A{v,i) = 16irv a s (-i)D(-i) 



«,(q 2 ) 



D-\c?) = c? + m\c?) 

12tt 



(33- 2^/) log 



q 2 +4m 2 (q 2 ) 
A 2 



(13a) 



(136) 



where the running gluon mass is given by 



2/ 2\ 2 

m (q ) = m G 



-, q 2 +4mn 

lOg A 2 ° 

1 4m n 

log^ 1 



-12/11 



(13c) 



The fixed mass ra^ can be determined PI from the condition that the simple exchange of a pair of 
gluons between quarks is the zeroth-order approximation to soft pomeron exchange at t = 0. This 
requires that the integral 

/3 2 = 1| d 2 qMq 2 )D(q 2 )] 2 (14) 
be about 4 GeV -2 . With a choice of A = 200 MeV this gives m = 340 MeV. 

While, strictly speaking, Cornwall's analysis applies only to Euclidean q, the form of (13) suggests 
that it effectively gives the gluon a mass which at low momentum values is close to m . Thus, while 
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Figure 4: Probability distribution for fractional energy of jets in the Born approximation 

there is uncertainty about how to calculate gluon emission, it seems that the confinement effects 
that Cornwall's calculation reveals do not allow the emission of soft gluons. So, while the effect of 
gluon emission tends to cancel that of gluon exchange, in the case of soft exchange there is no such 
cancellation. For this reason, we calculate the contribution from exchange for which q 2 < m 2 ,, and 
suggest that this provides a lower bound to the correction to the cross section from colour-octet 
exchange and emission. While we cannot be certain that this is a meaningful approach, we observe 
that the calculation is a very different one from that of the corresponding correction to the e + e~ total 
cross section. In the WW calculation there are W propagators in addition to quark propagators, which 
damp the contribution to the integration (3b) from large values of a and /3, so that the momentum 
transfer A 2 carried by the exchanged gluon is confined to spacelike values, as is seen in (4a). In the 
e + e~ case there is no such damping, and hence there are additional contributions from timelike A 2 , 
which we cannot calculate because the Cornwall form (13) of the gluon-exchnage amplitude is valid 
only for spacelike A 2 . 

We calculate the sum of the amplitudes of figure lb with the soft-gluon exchange attached to the 
quarks in all of the four possible ways*. We square this amplitude: colour considerations forbid 
interference with the Born term of figure 1, and of course in the squared amplitude there appears 
a colour factor 2/9 because the two gluons together must form a singlet configuration. We use the 
form (13a) for A in (9). We simply set fJ-i = fJ-2 = 0, because the nonperturbative quark propagator 
corresponding to the Cornwall gluon propagator is not available. The output for the ratio 1Z, defined 
in (11), is symmetric under x — > (1 — x) or y — > (1 — y) (or both). Figure 3 shows the contribution to 
1Z from the exchange of nonperturbative gluons with q 2 < 0.1 GeV 2 . The energy is y/s = 175 GeV 
and the plots are for Mi = M, with Mi = M and M ± V. The rather violent dependence on x and y 
is striking. 

The interaction has a significant effect on the energy distribution of the jets, though the overall effect 
on the integrated cross section is not so large. To convert the plots of figure 3 to energy distributions, 
note that the energy of the jet p\ is \X^fs, with 



* We do not consider diagrams where the gluon is exchanged simply between a pair of quarks associated 
with the same W, since this is just the same as the familiar correction to R e + e - and so is known to 



X = 



x(l - 2Ai + AiA 2 ) + Xi- AiA 2 
l-AiA 2 



(15) 



be small. 
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There are similar equations for the fractional energies of the other three jets. In the Born approxi- 
mation, the probability distributions of the two pairs of jets are uncor related, and each is symmetric 
under X — > (1 — X). Figure 4 shows the output, slightly smoothed, of a Monte Carlo calculation^] 
of the probability distribution P{X) at energy y/s = 172 GeV. We use this distribution to weight the 
output 1Z defined in (11) and we then average it over x and y. The result is shown in figure 5. 
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Figure 5: (TV) for y/s = 175 GeV, M 1 = M, plotted against M 2 . 

Conclusions 

We have found that colour transparency causes colour-singlet exchange to have very little effect on 
the cross section and spectra for the process e + e~ — > WW — > qqqq. We have calculated this by 
modelling colour-singlet in terms of a phenomenological soft-pomeron exchange, in which the soft 
pomeron couples to quarks with a form factor that goes to when the quarks go far off shell. We have 
shown elsewhere^] that when soft-pomeron exchange is modelled in terms of the exchange of a pair 
of nonperturbative gluons, this form factor arises from taking account of the various attachments of 
the gluons to the separate quarks of a colour-singlet system, that is it arises from colour-transparency 
effects. 

For colour-octet exchange the situation is very different. We calculate the contribution to this from 
soft-nonperturbative-gluon exchange, which we argue to give a lower bound to the total colour-octet 
exchange. Figure 5 shows that this will make a small but noticeable change in the total rate for the 
qqqq final state. In some regions of phase space the local effect can be significant, as can be seen from 
figure 3, and the energy distribution of jets will differ from that of the Born term alone. However 
there will be little impact on the mass determination due to the near symmetry of figure 5, although 
some increase in the width can be anticipated. 
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